Highly entangled quantum systems in 3 + 1 dimensions 
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Many systems exhibit boundary law scaling for entanglement entropy in more than one spatial 
dimension. Here I describe three systems in 3 + 1 dimensions that violate the boundary law for 
entanglement entropy. The first is free Weyl fermions in a magnetic field, the second is a holographic 
strong coupling generalization of the Weyl fermion system, and the third is a strong topological 
insulator in the presence of dislocations. These systems are unified by the presence of a low energy 
description that includes many gapless 1 + 1 dimensional modes. I conclude with some comments on 
the search for highly entangled states of quantum matter and some potential experimental signatures. 



Introduction. — The structure of entanglement in many 
body systems provides a useful characterization of these 
systems. It underlies the functioning of tensor network 
approaches to many body systems [lrl3j ■ It is also related 
to experimentally accessible quantities like current fluc- 
tuations Entanglement entropy, denned as the von 
Neumann entropy of the reduced density matrix of a cho- 
sen subsystem, provides a convenient characterization of 
entanglement. The emerging picture is that most many 
body systems live in a small subset of Hilbert space char- 
acterized by a geometric boundary law for entanglement 
entropy. Holographic duality has also provided insight 
into the deeply geometrical nature of entanglement [a- 
9] . I am interested in highly entangled states of quantum 
matter that violate the typical boundary law for the en- 
tanglement entropy. 

Many systems, gapless or not, satisfy a boundary law 
for entanglement entropy in higher dimensions (loj . This 
boundary law means that the entanglement entropy of a 
region of linear size L in d spatial dimensions scales as 
h d ~ x . There are exceptions. Systems described by con- 
formal field theory in one spatial dimension are known to 
violate the boundary in a universal way ll|, |12[ . The sim- 
plest exception in more than one dimension is provided 
by free fermions which violate the boundary with an ex- 
tra logarithmic correction [I3-[l3|- I have argued that 
Fermi liquids also violate the boundary law in a univer- 
sal way [1 81 ] - This argument is based on earlier intuition 
ascribing the boundary law violation of free fermions to 
the gapless one dimensional nature of excitations on the 
Fermi surface fl9j|. 

Here I use similar intuition to argue for boundary law 
violations in a variety of systems in higher dimensions. I 
will compute the entanglement entropy for three systems: 
a Weyl fermion in a magnetic field, a holographic system 
of many interacting Weyl fermions in a magnetic field 
[20| . and some stro ng top ological insulators with finite 



dislocation density [2lM27| . These systems all share a 



common feature which unifies the discussion, namely the 
appearance of a large number of gapless one dimensional 
modes at low energies. It is these one dimensional modes 
which are responsible for the violations of the boundary 
law. 



This paper is organized as follows. I begin by com- 
puting the entanglement entropy of a single free Weyl 
fermion in a background magnetic field. Next, I embed 
this Weyl fermion into a strongly coupled quantum field 
theory with a holographic dual and compute the entan- 
glement entropy for that system. I also describe a related 
setup in strong topological insulators. I conclude with a 
discussion of these results in the context of searches for 
highly entangled states of quantum matter. 

Weyl fermion in a magnetic field. — Consider a sin- 
gle Weyl fermion tp charged under a gauge field A with 
charge q in 3 + 1 dimensions. The equation of mo- 
tion for this fermion is j^D^ip — with j^ip = —i/j 
where = — iqA^ is the covariant derivative. Let 
there be a finite magnetic field, say in the z direction: 
F12 = d\A2 — 82 A\ = B. The magnetic field defines a 
length scale called the magnetic length l 2 B = 1/B (the 
units are made up by the flux quantum). On length 
scales much less than is the theory looks like a 3 + 1 di- 
mensional conformal field theory. On length scales much 
bigger than £b the theory becomes effectively 1 + 1 di- 
mensional. Indeed, the Weyl fermion is special because 
it possesses zero modes that avoid being gapped by the 
magnetic field. 

I take the gamma matrices to satisfy {7**, 7"} = 2n fil ' 
with rf u mostly minus. The chiral gamma matrix is de- 
fined to be 7 5 = Z7°7 1 7 2 7 3 and I work in the chiral basis 
where Dirac spinors decompose as ip T 
with 



7 5 = 



-1 2 
1 2 



The Weyl equation for a left handed spinor is 
(id t - ufD^L = 0, 



(1) 



(2) 



with a 1 the usual Pauli matrices. The vector potential in 
Landau gauge is A y = Bx for a constant magnetic field 
B in the z direction. Most solutions of the Weyl equa- 
tion in a finite magnetic field have a gap coming from the 
cyclotron motion, but there are also zero mode solutions. 
These solutions may be heuristically understood as aris- 
ing from a balance between the Zeeman energy and the 
orbital cyclotron energy. 
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Zero mode solutions may be found by putting dtipL — 
d z tpL = to obtain 

a x d x ^ L + <J y (d y - iqBx)^ L = 0. (3) 

Landau gauge maintains translation invariance in the 
y direction, so I try a solution of the form iPl( x i V) = 
ipL(x)e lky . The Weyl equation reduces to 



d x i(>L = -<J z {qBx - k)ip L 



(4) 



with solution 



V^)=exp(-^(z-^) o'jMO). (5) 

In order for this solution to be normalizable we must 
have ct z Vl(0) = ?/>z(0) (assuming qB > 0) leaving only 
one degree of freedom. The spacing of k is determined by 
the length of the system in the y direction to be Afc = y£. 

Lly 

We have one zero for each value of fc such that ip^ (x) sits 
inside the system in the x direction. The degeneracy g of 
zero modes is thus g — = g£_gia _ More generally, 

these zero modes and their degeneracy are protected by 
an index theorem relating the number of zero modes to 
the magnetic flux penetrating the system: iV zcro mo des — 
£jF 12 dxdy. 

So far we have ignored the z direction, but these 
zero modes actually disperse in the z direction. Assum- 
ing a more general solution of the form tpL(x,y, z,t) — 
e ip^z+t Py y-iEt q p L ^ full Weyl equation becomes 

E^l - p z v z TpL - i{o x d x + a y (d y - iqBx))ij L = 0. (6) 

The second half of this equation is solved with the same 
zero mode profile as above. The first half reduces to the 
equation E = p z using the fact that a z ^L = i>h follow- 
ing from the normalization condition. Thus each zero 
mode is actually relativistic chiral fermion in one spatial 
dimension. The low energy physics is controlled entirely 
by these zero modes as all other modes are gapped by 
the cyclotron motion. 

Using the one dimensional structure I can compute the 
entanglement entropy of the Weyl fermion. Consider a 
box of linear size L. The entanglement entropy Sl is de- 
fined as the von Neumann entropy of the reduced density 
matrix corresponding to the box: Sl — — Tr(p^ hip^). 
For one dimensional conformal field theories the entan- 
glement entropy is known to have the form 



Sr = 



CL 



(7) 



where cl and cr are the left and right central charges 
and e is an ultraviolet cutoff [Hj]. Weyl fermions in a 
magnetic field may be described by a large number of 
one dimensional gapless modes, and these modes are each 
equivalent to a chiral 1 + 1 dimensional conformal field 



theory, the dimensions being z and t. Each chiral fermion 
mode has cr, = 1 and cr = and hence contributes 
(1/6) In L to the entanglement entropy. For a cube of side 
length L aligned with the z direction we have qBL 2 / (2%) 
zero modes for a total entanglement entropy 



Sr 



( qBL 2 
V 2tt 



6 V e 



(8) 



This formula may be checked using the generalization of 
one dimensional entanglement entropy to finite temper- 
ature 



6 \7re p J 



(9) 



The thermal entropy of these zero modes in a cube of size 
L is thus 



S = 



qBL 2 \ irLT 



2tt J 6 



(10) 



which agrees with the direct thermodynamic calculation. 

Before moving on, let me note that a single charged 
Weyl fermion does not give a consistent quantum theory. 
This is due to the presence of a gauge anomaly propor- 
tional to Tr(Q 3 ) where Q is the charge matrix. There is 
also a gravitational anomaly proportional to Tr(Q). Both 
of these anomalies must vanish for a completely well de- 
fined chiral gauge theory, but this can be accomplished by 
adding Weyl fermions with compensating charges. The 
boundary law violating behavior remains, and thus there 
are consistent configurations of Weyl fermions that vio- 
late the boundary law for entanglement entropy. 

Holographic generalization. — I have computed the en- 
tanglement entropy for a single free Weyl fermion and 
found a term that violates the boundary law for entangle- 
ment entropy. A useful choice for incorporating interac- 
tions is J\f = 4 SU(A r ) Yang-Mills theory which includes 
4iV 2 Weyl fermions as part of the field content. These 
fermions sit in the adjoint of the non-Abelian gauge group 
SU(A r ), while the magnetic field B corresponds to a 
weakly gauged U(l) subgroup of the R-symmetry. In 
zero magnetic field this theory is conformal at all values 
of the t'Hooft coupling A = gy M N, but it is particularly 
amenable to study at strong coupling because of holo- 
graphic duality. This duality relates the N — 4 theory 
to a theory of quantum gravity, IIB string theory, in an 
asymptotically five dimensional anti-de-Sitter spacetime 
(AdSs). The limit A — > oo and N — » oo in the field the- 
ory gives classical supergravity in anti-de-Sitter space on 
the gravity side. 

In this strong coupling limit, configurations of the su- 
per Yang-Mills theory have an emergent geometric inter- 
pretation in terms of classical gravitational field config- 
urations. The ground state of the field theory is dual to 
pure anti-de-Sitter space, and the field theory at finite 
temperature is accessed via a bulk black hole. The field 
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theory in a background magnetic field at zero temper- 
ature is obtained from a magnetically charged extremal 
black hole in the bulk. Given the bulk geometric config- 
uration, the leading large N contribution to the entan- 
glement entropy can be determined holographically by 
computing the area of certain minimal surfaces in the 
bulk @,|J. 

Consider extremal magnetic brane solutions in 
Einstein-Maxwell theory with negative cosmological con- 
stant in five dimensions [20] . These solutions interpolate 
between an asymptotically AdSs region and a near hori- 
zon AdS3 x T 2 region (assuming the xy plane is com- 
pactified). The asymptotic AdSs region corresponds to 
the unperturbed Af = 4 theory at high energies. The 
near horizon region appears as a result of turning on a 
magnetic field in the gauge theory. The radial evolotion 
represents a renormalization group flow from a 3 + f di- 
mensional conformal field theory at high energies to an 
effectively 1 + 1 dimensional conformal field theory at 
low energies. This is qualitatively similar to the physics 
of free Weyl fermions, and even at strong coupling the 
cross-over scale is determined by the magnetic length. 
At zero temperature the metric may be written in the 
form 



ds 2 = -U(r) dt 2 + 



dr 2 



+ U{r) dz l + e 2V {dx z +dy z ), 



(11) 

with r the radial coordinate (r — ► oo is the boundary) 
and z the direction of the magnetic field on the boundary 
[20j |. I use bulk units with the AdS radius set to one. 
In addition to the metric, the gauge field has a profile 
given by F — B dx A dy. The asymptotic AdSs region 
is described by U = e 2V = r 2 while the near horizon 
AdS 3 x T 2 region corresponds to U = 3r 2 and e 2V = B/3. 
Notice that in the near horizon region the xy plane has 
decoupled from the radial coordinate and has fixed size 
given by the magnetic length. 

The entanglement entropy of a region in the dual field 
theory is determined by the area of the minimal surface 
in the bulk that terminates on the boundary of the re- 
gion in the field theory. The entanglement entropy is just 

(5) 

this minimal area divided by AG N . I will focus on the 
entanglement entropy of a rectangular region in bound- 
ary theory of size L x L x L z . Assuming L 3> L z gives 
approximate translation invariance in the xy plane. The 
minimal surface calculation reduces to a two dimensional 
problem involving only the variables z and r. The zero 
temperature geometry is only known numerically, and the 
minimal surface calculation can also only be done numer- 
ically. However, the important physics can be extracted 
without the numerical details. For cubic regions with all 
dimensions less than the magnetic length, the minimal 
surface only probes the AdSs region and gives the usual 
ultraviolet divergent boundary law for entanglement en- 
tropy. 



For boundary regions of linear size much larger than 
the magnetic length, the minimal surface passes right 
through the asymptotic AdSs region towards the near 
horizon region. Once in the near horizon region, the 
x and y directions freeze out, and the minimal surface 
behaves exactly as in AdS3. In particular, I find the 
characteristic \n(L z /£) dependence familiar from 1 + 1 
dimensional conformal field theory with the magnetic 
length providing the cutoff. The entanglement entropy 
thus consists of two pieces, a non-universal boundary law 
contribution from the asymptotically AdSs region and a 
universal low energy piece Sl ~ N 2 BL 2 In (L z /£b)- The 
appearance of the magnetic field can be understood be- 
cause the effective 1 + 1 dimensional central charge is 
related to l/G$ which is enhanced relative to 1/G^ by 
a factor of BL 2 from the freeze out of the xy plane. This 
strong coupling version of the free Weyl fermion system 
thus also violates the boundary law for entanglement en- 
tropy at low energies. 

Strong topological insulators. — In both the cases con- 
sidered above, the appearance of many gapless one di- 
mensional modes was responsible for the highly entan- 
gled nature of the quantum state. This intuition can be 
applied to more experimentally relevant systems known 
as strong topological insulators. These systems are time 
reversal invariant electronic band insulators that are not 
smoothly connected to trivial band insulators. In partic- 
ular, they possess interesting topological structure that 
gives rise to protected edge modes. These edge modes are 
robust so long as time reversal invariance is preserved 
26|. Topological insulators in three spatial di- 



21, 23 



mensions have gapless surface states living in two spatial 
dimensions, but these modes do not lead to a violation 
of the boundary law for entanglement entropy. Similarly, 
the bulk of a topological insulator is gapped in a perfect 
crystal and certainly satisfies a boundary law for entan- 
glement entropy. 

However, experimentally realized topological insula- 
tors are not perfect crystals, they possess topological de- 
fects including dislocations in the crystalline bulk. Re- 
markably, for certain kinds of topological insulators and 
dislocation types, the dislocations have been shown to 
support gapless fermionic modes [28|. These effectively 
one dimensional modes make the dislocations into gapless 
quantum wires threading the otherwise gapped bulk. The 
one dimensional modes in the quantum wires are analo- 
gous to the Weyl zero modes considered above, with the 
dislocations playing the role of magnetic field lines. In 
the presence of a finite density of dislocations supporting 
gapless modes, the bulk of a strong topological insulator 
violates the boundary law for entanglement entropy. 

To estimate the size of the violation, consider the artifi- 
cial situation of a dilute array of topologically non-trivial 
dislocations all aligned. Let these dislocations have an 
areal density p (a typical value of p might be 10 12 m~ 2 
28]). A region in the bulk of size L X L X L z , with 
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the z axis chosen parallel to the dislocations, effectively 
contains pL 2 gapless one dimensional fermionic modes. 
These modes should each contribute roughly In (L z /e) to 
the entanglement entropy. The boundary law violating 
component of the entanglement entropy is thus of order 
Sl ~ pL 2 In (L z /e). This estimate is crude, but it should 
suffice for a reasonably uniform and collimated set of dis- 
locations. Note that despite the enhanced L dependence 
relative to the usual boundary law, this term may be 
much smaller than the boundary law term for experimen- 
tally accessible system sizes and dislocation densities. I 
also wish to emphasize that this is a statement about 
the zero temperature quantum state. The helical modes 
modes are protected from elastic scattering (such scat- 
tering might otherwise localize a one dimensional gapless 
mode), but at finite temperature or in the presence of 
inelastic processes, the boundary law violating behavior 
will be disrupted. 

Discussion. — I have described three systems in 3 + f 
dimensions that violate the boundary law for entangle- 
ment entropy. The mechanism for these violations is the 
emergence of a low energy description in terms of many 
gapless one dimensional modes. A completely analogous 
picture has been argued to hold for Fermi liquids in any 
dimension, with the exact coefficient of the boundary law 
violating term determined by the geometry of the inter- 
acting Fermi surface. The generic picture of modes prop- 
agating along one dimensional objects, be they magnetic 
field lines or topological defects, permits a unified under- 
standing of the weakly coupled systems I considered. In 
the strong coupling limit, holographic duality provides 
a geometric realization of a similar emergent 1 + 1 di- 
mensional description in the magnetic field direction and 
time. 

These systems are interesting in part because they ex- 
pand our understanding of the limits of the boundary 
law for entanglement entropy. Currently, violations of 
the boundary law in reasonable (without large ground 
state degeneracy) translation invariant systems can be 
traced to a one dimensional picture. The one dimen- 
sional physics in these systems suggests the possibility of 
enhanced fluctuations associated with violations of the 
boundary law @, H, HI] . A different test might come 
from thermodynamic measurements. The thermal signa- 
ture of boundary law violating entanglement entropy in 
these effectively relativistic one dimensional systems is a 
contribution to heat capacity linear in T . In the case of 
topological insulators, this term will be quite small, but 
it does have a different temperature dependence than the 
contributions from surface states, gapped electrons, and 
phonons. 

There are many directions for future work. In the di- 
rection of classification, it is an interesting open question 
whether all violations of the boundary law in reasonable 
translation invariant systems can be traced to one dimen- 
sional physics. Are there examples of reasonable quan- 



tum systems even more highly entangled than these ef- 
fectively one dimensional systems? On the experimental 
front, it has proved difficult to find definite signatures 
of entanglement in many body systems in higher dimen- 
sions. The situation is better in one dimension where 
the entanglement entropy can be directly measured using 
current fluctuations Q ■ It would be interesting to explore 
in detail the experimental consequences of entanglement 
in higher dimensional systems that are nevertheless con- 
trolled by one dimensional physics. 
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